Propagation of exciton pulses in semiconductors 
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Using a toy model, we examine the propagation of excitons in CU2O, which form localized pulses 
under certain experimental conditions. The formation of these waves is attributed to the effect of 
dispersion, non-linearity and the coupling of the excitons to phonons, which acts as a dissipative 
mechanism. 
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I. INTRODUCTION 

The transport properties of excitons [1] in CU2O have 
been studied extensively in recent years. In one series of 
experiments [2-6] , a pulse of excitons was created at one 
end of a crystal of CU2O and detected at the other fol- 
lowing propagation inside the sample. In another series 
of experiments [7-9], excitons were created and imaged 
as they expanded inside a crystal subjected to external 
stress. (The stress effectively created a confining poten- 
tial for the excitons.) 

Two noteworthy conclusions emerged from these ex- 
periments. First, drag in the propagation of excitons due 
to their scattering with (longitudinal) acoustic phonons 
was found to decrease with decreasing lattice tempera- 
ture. The second observation [3-5] followed from experi- 
ments with two laser pulses separated by a time delay of 
~ 100 ns. The first pulse created an asymmetric wave of 
excitons; the excitons created by the second pulse trav- 
elled in the trail of the first. This second wave was ob- 
served to be symmetric and to have the characteristics 
of a solitary- wave, specifically, those of a "bright soliton" 
(i.e., a local elevation in the density.) 

Here, we will argue that the waves created by the sec- 
ond pulse are likely to resemble classical solitary waves. 
These waves are, however, essentially different from soli- 
tary waves in, e.g., water (as first observed by J. Scott 
Russell [10]) due to the presence of dissipation caused by 
the scattering of excitons with phonons [7-9]. We will 
suggest that the effects of nonlinearity and dispersion, 
which give rise to solitons, must be combined with dissi- 
pation in order to create generic conditions appropriate 
for the propagation of localized pulses of excitons in the 
crystal. All three mechanisms play an important role. 

The studies of Refs. [11,12] attributed the observed 
localized waves to the formation of a Bose-Einstein con- 
densate of excitons [13]. In Ref. [11], second-order per- 
turbation theory was used to describe the exciton-exciton 
interaction [11], which gave rise to effective attraction 
between the condensed excitons. In Ref. [12], it was 
shown that an exciton condensate coupled to acoustic 
phonons can give rise to bright solitons. However, as 



pointed out in Ref. [14], solitary waves in a Bose-Einstein 
condensed gas have a characteristic length scale deter- 
mined by its healing length £ given by Ti /2m£ 2 = hUq, 
or £ = (87rna) -1 ' 2 , where m is the exciton mass, n is the 
exciton condensate density, and Uq = AttTi a/m is the 
exciton-exciton scattering matrix element, with a being 
the scattering length. For typical liquid helium temper- 
atures of order 2 K, the critical density for the Bose- 
Einstein condensation of excitons is of order 10 17 cm~ 3 . 
The scattering length for exciton-exciton elastic collisions 
in Cu 2 has been calculated to be a « 2clb [15], where 
as is the exciton Bohr radius. The Bohr radius of exci- 
tons in CU2O is rs 5 A [16] so that £ ~ 0.01 /im, which is 
about four orders of magnitude smaller than the width of 
the narrowest pulses l p ~ 1 mm. (The crystal is typically 
a few millimeters long.) 

Finally, a model which has been proposed in connec- 
tion with the ballistic propagation of excitons attributes 
the observed phenomenon to a classical effect due to a 
"phonon wind" [17]. Many of the observations have been 
explained in the context of this model and our study is 
consistent with the claim that the excitons are described 
by classical equations. On the other hand, the two stud- 
ies examine effects which are essentially decoupled. 

More specifically, while both models include dissipa- 
tion, the phonon-wind model focuses on the spatial and 
temporal effects of laser excitation, whereas our model 
concentrates on the effects of non-linearity and disper- 
sion. In the present study we make the implicit simpli- 
fied assumption of a delta- like form of laser excitation in 
space and time. In such a case, the phonon-wind model 
essentially reduces to a diffusion-like equation, neglecting 
the effects of non-linearity and dispersion. Of course, one 
should worry about more complicated forms of the laser 
excitation, however the problem we have studied is quite 
complicated and we have chosen to overlook these ef- 
fects in order to demonstrate clearly the combined effect 
of dissipation, non-linearity and dispersion; we believe 
all three effects are very crucial in these experiments. 
Finally, from an experimental point of view, our study 
would be mostly suitable for short (in space and time) 
laser pulses. 



II. HYDRODYNAMIC DESCRIPTION 

Following Rcf. [18], we start with the hydrodynamic 
equations describing the system, 
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where v is the velocity field, p is the pressure, and r is 
the exciton-acoustic phonon scattering time. The neglect 
of viscosity in these equations can be justified by looking 
at the characteristic scales: The viscous force nW 2 v is of 
order nv/l 2 ~ nT s c 3 /l 2 , where n is the viscosity, t s is the 
scattering time for elastic exciton-exciton collisions, and 
c 2 = dp/dn is the adiabatic sound velocity. Therefore, 
\n\7 2 v/Vp\ ~ cT s /lp. Typical values of c ~ 5 x 10 5 cm/s, 
l p ~ 1 mm, and r s = (nav)^ 1 with a = 8ira 2 imply 
that t s ~ 1 ns for n ~ 10 17 cm~ 3 , a ~ 10 A, and a 
temperature T ~ 1 K. Therefore, ct s /I p ~ 10~ 2 <C 1. 

Linearizing the above equations and considering small 
variations in the density and the pressure, n = no + <5n 
and p = po + Sp (where no and po & r e the "background" 
density and pressure of excitons created by the first laser 
pulse), one can derive the equation: 
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dt 2 t dt 

This equation leads to the dispersion relation [18] 

uj 2 - c 2 k 2 + iiu/T = 0, (3) 

which interpolates between two limits. For kcr <C 1, the 
motion of the exciton pulse is diffusive, and 

d t Sn = c 2 T\7 2 5n, (4) 

with u) = —ic?Tk 2 . In the opposite limit kcr ^> 1, 

d tt Sn - c 2 \7 2 5n = 0, (5) 

with u) — ck, corresponding to sound propagation. 

III. NON-LINEARITY VERSUS DISPERSION 

The analysis so far applies to small variations around 
equilibrium. For larger excitations it is also necessary 
to account for the effects of non-linearity and dispersion. 
For simplicity, we start with a (one-dimensional) wave 
propagating with velocity c along the positive x axis, 



dtSn + cd x Sn = 0. 
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Any function of the form Sn — 5n{x — ct) is a solution 
of Eq. (6). Allowing for dispersion and non-linearity, this 
becomes 
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Here a and /? are small (positive) dimcnsionless quanti- 
ties, which can be calculated from microscopic theories. 
It is not our goal to calculate them here, and we treat 
them as phenomcnological parameters. Also t p is the 
time for the pulse to travel a length l p (i.e., l p = cr p .) 
Typically, r p ~ 100 ns for relatively narrow pulses. Equa- 
tion (7) leads to the dispersion relation 
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which provides a reasonable description of non-linearity 
and dispersion. For k — > 0, the sound velocity is c(l + 
aSn/no), and the local sound velocity is greater in regions 
of higher density. For a = 0, uo/k = c[l — f3(kl p ) 2 ], which 
implies that longer wavelengths travel faster. Equation 
(7) is of standard Korteweg-de Vries form, and the com- 
petition between the last two terms gives rise to solitonic 
solutions [19]. 

We seek solutions of Eq. (7) which preserve their shape, 
i.e., which have the form 5n(x — ui), and which satisfy 
the boundary condition that Sn — > as \z\ — > 00. Such 
solutions have the form 
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with u = c(l + aC/3). For fixed a and /?, there is thus a 
family of solutions characterized by their maximum am- 
plitude, C. For aC ~ 12/3 the size of the travelling wave 
is of order l p . This result reflects the expected competi- 
tion between non-linearity and dispersion that gives rise 
to this localized excitation. In addition, the velocity of 
this pulse is of order c, which can approach the speed 
of sound ci of CU2O. For example, for a classical gas 
p = nksT ', and thus c 2 = ksT '/m, with m = 3m e being 
the total exciton mass density (m e is the electron mass.) 
For T w 2 K, c w 3.2 x 10 5 cm/s, in good agreement with 
experiment. 



IV. SCATTERING OF EXCITONS WITH THE 
CRYSTAL 

Equation (7) docs not include the effects of the scat- 
tering of excitons with the crystal. As noted above, this 
mechanism leads to the diffusive motion of excitons. It 
is thus worth considering what is known about exciton- 
phonon interactions in CU2O. 

In fiefs. [7,8], paraexcitons under low excitation condi- 
tions were observed to be highly mobile, with a diffusion 
constant D ~ 10 3 cm 2 /s. For D = c 2 t with c 2 = ksT/m, 
t is seen to be w 20 ns. For typical thermal velocities 
of order 3 x 10 5 cm/s, the mean-free path of excitons is 
w 50 /im. This corresponds to w 10 5 lattice sites and 



is thus quite large. In addition, D showed a rapid in- 
crease for temperatures below w 6 K. These observa- 
tions were explained in Ref. [8] using a simple classi- 
cal model which showed that excitons with wavevectors 
smaller than ki — mci/h, where q is the longitudinal 
speed of sound in CU2O (cj « 4.5 x 10 5 cm/s), can- 
not emit acoustic phonons and thus cannot scatter with 
the lattice. For exciton temperatures below 7], where 
^bTi ~ mcf, exciton-acoustic phonon scattering is sup- 
pressed. This freeze-out temperature is TJ « 4 K for lon- 
gitudinal acoustic phonons in Cu 2 0. 

In the experiments of Rcfs. [2-5] with an excitation of 
w 1.2 x 10 6 W/cm 2 , the temperature varied between 1.9 
K and 4.1 K. For T w 2.5 K, which is close to the freeze- 
out temperature Ti, the propagation of the exciton pulse 
showed a transition from diffusive to ballistic motion. In 
Rcf. [4] for conditions under which the exciton propa- 
gation was diffusive, the diffusion constant was reported 
to be even higher, Dw2x 10 3 cm 2 /s, corresponding to 
t sw 40 ns. 

It is worth pointing out that in the experiments where 
the sample was illuminated with two pulses [3-5], these 
were delayed by ~ 100 ns. This time delay is compara- 
ble to t p . It is thus reasonable to expect that the first 
pulse diffuses inside the crystal and creates a background 
density of excitons through which the second pulse must 
pass. 



V. A MORE REALISTIC PROBLEM 

In view of these considerations, we have included 
an additional term in Eq. (7) to describe the effects of 
exciton-phonon scattering: 

dSn , Sn s dSn n L d 3 5n l„ d 2 5n , , 

■flr + c(1 + ^-to + ^ft?r = ^-^' (10) 

with 7 being a positive parameter. Here t p /"/ can be 
identified as the exciton-acoustic phonon scattering time. 
Since t p ~ 100 ns and r is ~ 30 ns, 7 is expected to be 
greater than, but on the order of unity. The term on 
the right side of Eq. (10) gives an imaginary part to the 
dispersion relation, modifying Eq. (8) to 
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With such a dispersion relation the wave is damped. For 
example for k ~ l/l p , Sn ~ e~ 7 */ Tp . 



velocity of order c. When 7 is nonzero, the width of the 
pulse increases with time and its shape changes. As 7 
is increased, this spreading velocity increases and even- 
tually equals c at some 7 max - For 7 <C 7max, the pulse 
pulse spreads slowly relative to its propagation velocity 
of « c. 

To determine 7 maX j it is convenient to use the notation 
(f) = Su/uq and the quantities 



(X) = 
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and (X 2 ) = 
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and the associated width AX = ^/{X 2 ) - (X) 2 . The 
time derivative d t (X) can be regarded as the velocity 
of propagation of the pulse; d t AX describes the rate of 
increase of the width. The desired average values follow 
immediately from appropriate integrals of Eq. (10). We 
find that 
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For 7 = and <j> given by Eq. (9), Eq. (13) gives the ex- 
pected result that dt(X) = c(l + aC/3), i.e., the velocity 
of propagation u given above. Similarly, we find that 
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The a-dependent term in this equation evidently vanishes 
for any symmetric pulse shape. For such pulses we see 
that 



d f AX = 



AX' 



(15) 



For pulses of typical width of order l p , this suggests that 
7max is of order 1. The characteristic scattering time 
of « Tp, which establishes the crossover from the mo- 
tion of a localized pulse to diffusive motion, is thus of 
the same order as the characteristic scattering time r for 
exciton-phonon collisions that has been determined ex- 
perimentally as described above. As we discuss below, 
the approximately linear growth of AX 2 (t) indicated by 
Eq. (14) provides a definite test for the validity of our 
description. 



VII. THE BENEFITS OF DISSIPATION 



VI. EFFECT OF DISSIPATION ON THE PULSE 
PROPAGATION 

Equation (10) interpolates between two extreme modes 
of propagation. For 7 = the solution is a travelling 
wave, which propagates without change in shape with a 



While dissipation can obscure the solitonic propaga- 
tion, it also plays a far more positive role. In its ab- 
sence, an initially localized pulse will often break up into 
a discrete number of solitons (plus small amplitude, non- 
solitonic motion.) This is illustrated in Fig. 1(a) where 
we show the evolution of an initial pulse 
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with A = 3, which leads to the formation of three solitons. 
(The calculation was performed for a = 0.24, j3 = 0.02, 
and 7 = 0.) Once formed, the distance between the indi- 
vidual solitons grows linearly with t. The generic break- 
up of localized peaks into many solitons is in marked con- 
trast to the experiments considered here, which see only 
single, localized peaks. Dissipation of sufficient strength 
can provide an effective mechanism for preventing the 
formation of unwanted local maxima. This is illustrated 
in Figs. l(b)-(d) obtained numerically for 7 = 0.01, 0.03, 
and 0.2, respectively. Once new local extrema have been 
formed, dissipation will tend to reduce the local maxi- 
mum (where the curvature is negative) and enhance the 
local minimum (where the curvature is positive.) For 
all values of 7, the local maximum and minimum slowly 
merge and ultimately disappear. Given the linear growth 
of the separation between local maxima and the much 
longer time scale required for dissipation to eliminate 
them, it is clear that 7 should exceed some minimum 
value, 7 m in, chosen to ensure that secondary maxima 
never form. The value of 7 m i n is set by a complicated 
interplay between the initial pulse shape and the param- 
eters appearing in Eq. (10). For the present case of A = 3, 
numerical simulations indicate that 7 m i n = 0.024. The 
initial pulse width considered here is larger than that ap- 
propriate to the experimental initial conditions of Rcfs. 
[2-5]. A more realistic calculation would involve a value 
of A that is smaller than (but on the order of) unity. We 
have used here the value A = 3 for convenience. In addi- 
tion, numerical evidence indicates that 7 m i n grows slowly 
as the width A decreases. 

An interesting feature of our results is that 7 max and 
7 m in differ by a large factor (i.e., 40 in the present case.) 
This suggests that it is not necessary to fine tune the dis- 
sipative term in Eq. (10) in order to obtain localized pulse 
propagation. Indeed, the qualitative impression of soli- 
tary pulse propagation seems to be a particularly robust 
phenomenon, and its observation cannot be regarded as 
a demonstration of true solitonic propagation. 



VIII. TWO EXPERIMENTAL CHALLENGES 

The model we have presented can be checked eas- 
ily since it makes the following two definite predictions: 
1) We have noted that, in the presence of dissipation, 
AX 2 (t) grows linearly with t. Since t is the time for the 
pulse to propagate inside the crystal, AX 2 should also 
grow linearly with the size L of the crystal along the di- 
rection of pulse propagation (keeping all other conditions 
fixed.) 2) As we have seen, the propagation of excitons is 
sensitive to the value of 7, and thus to their coupling with 
phonons. If an external stress (of a few kbar) is applied, 
excitons can couple to transverse acoustic phonons, as 
discussed in Ref. [8] . The transverse speed of sound ct is 




FIG. 1. Solutions of Eq. (10). We have plotted Sn/n (y 
axis) as function of distance (x axis), measured in units of l p . 
The initial disturbance is given by Eq. (16) with A = 3. We 
have used the values a = 0.24 and /3 = 0.02. The figures show 
snapshots of the solution for 7 = (a), 0.01 (b), 0.03 (c), and 
0.2 (d), at times t/r p = 10.0, 25.0, and 50.0. 



lower than c; , and the corresponding freeze-out tempera- 
ture is consequently some ten times smaller than Tj, i.e., 
« 0.5 K [8]. If the same experiments were performed in 
the presence of stress, the term on the right of Eq. (10) 
would be enhanced and the motion would become more 
diffusive with increasing stress. (The conclusions of the 
other theories [11,12,17] are expected to be less sensi- 
tive to external stress.) These tests are likely to provide 
compelling evidence for or against our picture since the 



empirical value of 7, estimated above, is relatively close 
to 7 max . 



IX. SUMMARY 

We have presented a phenomenological theory at- 
tributing the formation of localized waves of excitons in 
CU2O under certain experimental conditions to the com- 
bined effects of dispersion, non-linearity and dissipation. 
Both analytical estimates and numerical simulations sug- 
gest that the characteristic time scale for the crossover 
between the regimes of localized propagation of the pulse 
and diffusive motion is comparable to the scattering time 
for exciton-phonon collisions. 
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